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Recent experimental developments of high-intensity, short-pulse XUV light sources are enhanc- 
ing our ability to study electron-electron correlations. We perform time-dependent calculations to 
investigate the so-called "sequential" regime (hco > 54.4 eV) in the two-photon double ionization of 
helium. We show that attosecond pulses allow to induce and probe angular and energy correla- 
tions of the emitted electrons. The final momentum distribution reveals regions dominated by the 
Wannier ridge break-up scenario and by post-collision interaction. 

PACS numbers: 32.80.Rm, 31.15.V-, 32.80.Fb, 42.50.Hz 



Understanding the role of electron correlation in 
atoms, molecules, and solids has been a central theme 
in physics and chemistry since the early days of quantum 
mechanics. Most of the focus has centered on the role 
of electron correlation in (quasi-) stationary states. Re- 
cent progress in the development of light sources provides 
unprecedented opportunities to expand our understand- 
ing of electron correlation to dynamical processes where 
external fields play a critical role. The availability of 
attosecond pulses, as generated from high-harmonic ra- 
diation [1-3], opens new avenues for time-domain studies 
of multi-electron dynamics. Using such pulses, it is pos- 
sible to not only observe, but also actively induce and 
control correlation effects. 

The simplest system where electron-electron interac- 
tion can be studied is the helium atom. Unraveling the 
intricacies of electron correlation in ultrashort and in- 
tense electromagnetic fields interacting with this simple 
atom is critical to our understanding of the same pro- 
cesses in more complex systems. Despite the computa- 
tional challenges, the dynamics of He under the influ- 
ence of external fields can still be accurately simulated in 
ab initio calculations, cf. [4]. The results of the present 
investigation provide evidence that the effects of elec- 
tron correlation can be surprisingly complex in situa- 
tions dominated by external ultrashort fields. This in 
turn has important consequences for attosecond studies 
in molecules, clusters, and solids. We show that it is pos- 
sible to disentangle the different processes occurring in 
such pulses by analyzing the final momentum distribu- 
tion of the ejected electrons. 

Double ionization of helium by single photon absorp- 
tion has long been the benchmark for our understanding 
of correlation effects in the three-body Coulomb problem 
[5-9] . The availability of intense light sources in the VU V 
and XUV region [10-12] has recently shifted attention 
from single-photon double ionization and intense-IR laser 
ionization by rescattering (see [13-15] and references 



therein) to multiphoton ionization. Restricting attention 
to only two-photon double ionization (TPDI) enables us 
to distinguish two spectral regions. The "nonsequential" 
or "direct" regime between 39.5 eV < hw < 54.4 eV has 
recently received considerable attention (see [16-21] and 
references therein). Energy-sharing between the elec- 
trons, and thus correlations, are a conditio sine qua non 
for double ionization to occur in this regime. By contrast, 
in "sequential" TPDI with hw > 54.4 eV [22-25], each 
photon has sufficient energy to ionize one electron within 
an independent-particle model and electron-electron in- 
teraction, while present, is not a necessary prerequisite. 

For an ultrashort pulse of attosecond duration the con- 
cept of "sequential interactions", valid for long pulses, 
becomes obsolete. Instead, the two-electron emission 
occurs almost simultaneously, and the strength of elec- 
tron correlation in the exit channel can be tuned by 
the pulse duration T p . This information is encoded 
in the final joint momentum distribution P(ki,k2) = 
P(Ei, £2, ^1, ^2), experimentally accessible in kinemat- 
ically complete COLTRIMS measurements [26]. 

In our current calculations we solve the time-dependent 
Schrodinger equation in its full dimensionality, including 
all inter-particle interactions. The laser field is linearly 
polarized and treated in dipole approximation. The du- 
ration T p is given by the FWHM of a sine-squared en- 
velope function for the electromagnetic field. The com- 
putational approach is based on a time-dependent close- 
coupling (TDCC) scheme where the angular variables are 
expanded in coupled spherical harmonics and the two ra- 
dial variables are discretized via a finite element discrete 
variable representation (FEDVR) . Temporal propagation 
is performed by the short iterative Lanczos (SIL) algo- 
rithm with adaptive time-step control. The asymptotic 
momentum distribution is obtained by projecting the 
wave packet onto products of Coulomb continuum states. 
Projection errors due to the replacement of the full three- 
body final state by independent-particle Coulomb wave 
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FIG. 1: TPDI electron spectra P(E U E 2 ) at hu = 70 eV 
for different pulse durations (FWHM): (a) T p = 150 as, (b) 
T p = 750as. The top shows the spectrum integrated over one 
energy, i.e., the one-electron energy spectrum P(E\) = P(E2). 



functions can be reduced to the one-percent level by de- 
laying the time of projection until the two electrons are 
sufficiently far apart from each other [19]. 
The joint energy probability distribution 



P(E U E 2 ) = JJ P(E u E 2 ,Sl u Sl 2 )dSl 1 dSl 2 



(1) 



reveals the breakdown of the sequential ionization pic- 
ture with decreasing pulse duration T p (Fig. 1). For 
long pulses, two distinct peaks signifying the emission 
of the "first" electron with energy E\ = huj — 1\ (with 
l\ = 24.6 eV the first ionization potential) and the "sec- 
ond" electron with E 2 = Huj — I 2 (22 = 54.4 eV) are clearly 
visible. 

For pulses of the order of one hundred attoseconds a 
dramatically different picture emerges: the two peaks 
merge into a single one located near the point of sym- 
metric energy sharing. It should be noted that this is 
not simply due to the Fourier broadening of the pulse. 
Instead, the close proximity in time of the two emis- 
sion events allows for energy exchange between the two 
outgoing electrons representing a clear departure from 
the independent-particle behavior [22, 25]. Differently 
stated, the time interval between the two ionization 
events is too short for the "remaining" electron to re- 
lax to a stationary ionic state. In the limit of ultrashort 
pulses the notion of a definite time ordering of emission 
processes loses its significance, as does the distinction be- 
tween "sequential" and "nonsequential" ionization. 

The attosecond-pulse induced dynamical electron cor- 
relation becomes more clearly visible in the joint angu- 
lar distribution P(0i2,#i) (Fig. 2), where 6\ is the polar 
emission angle of one electron with respect to the po- 
larization axis of the XUV pulse, 9\ 2 is the angle be- 
tween the two electrons, and the energies E\^E 2 are 
integrated over. Here and in the following we choose 
coplanar geometry with <pi = <p 2 = 0°. In the limit of 
"long" pulses (T p = 4.5 fs), the joint angular distribu- 
tion approaches the product of two independent Hertz 




FIG. 2: Conditional angular distributions P(0i2,0i=O°) of 
the ejected electrons for different pulse lengths at huj = 70 eV. 
The innermost (solid blue) line is for T p = 75as FWHM, with 
successive lines for T p = 150 as, 300 as, 750 as, and 4500 as 
FWHM. For better comparison the distributions are normal- 
ized to a maximum value of one. 
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FIG. 3: Conditional angular distributions P(#i2,#i=0°) for 
a 75 as (FWHM) pulse at /kj = 70eV (solid blue) and for 2fs 
(FWHM) pulses at 42 eV (dashed red) and 52 eV (dash-dotted 
green). The distribution for the ultrashort pulse strongly 
resembles the long-pulse distribution in the nonsequential 
regime (hco < 54.4 eV). 



dipoles, each of which signifies the independent interac- 
tion of one electron with one photon. Consequently, also 
the conditional angular distribution P(#i2,#i =0°) cor- 
responds to a Hertz dipole. With decreasing pulse du- 
ration, P(#i2,#i =0°) is strongly modified and develops 
a pronounced forward-backward asymmetry. The condi- 
tional probability for the second electron to be emitted 
in the same direction as the first is strongly suppressed. 
It is worth noting that this strong preference for back-to- 
back emission for T p < 150 as persists after integration 
over the electron energies. Nevertheless, approximately 
equal energy sharing dominates (cf. Fig. 1). Thus, the 
dominant break-up mode induced by an attosecond pulse 
corresponds to the "Wannier ridge" configuration [27]. 
The same break-up mode is observed in the nonsequen- 
tial TPDI regime (hco < 54.4 eV, cf. Fig. 3), where the 
electrons need to exchange energy to achieve double ion- 
ization. Thus, even in long pulses only electrons ionized 
within a short time of each other can be observed. 

It is now instructive to inquire into the origin of the 
strong angular correlations observed for short pulses. 
Three different sources can be distinguished: 

(i) Correlations in the helium ground state. Due to 
Coulomb repulsion, the electrons in the ground state are 
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FIG. 4: Conditional angular distributions, (a) For a duration 
(FWHM) of 75 as for different photon energies. From inside 
to outside: 70 eV, 91 eV, 140 eV, and 200 eV. The amount of 
asymmetry decreases with increasing pulse energy, (b) For 
different observation times after the 75 as FHWM pulse at 
hco = 70eV. Snapshots were taken (from the outermost to the 
innermost line) immediately at the end of the pulse and 50 as, 
150 as, 600 as, and 1000 as after the end of the pulse. 

not independent of each other. For ultrashort pulses, 
TPDI can thus be interpreted as a probe that maps out 
the initial-state correlations. 

(ii) Induced dipole polarization in the intermediate 
state. When the first electron leaves the core, its electric 
field induces polarization of the remaining ion, leading 
to an asymmetric probability distribution of the second 
electron. The second photon then probes the dynam- 
ics in this bound-free complex, such that TPDI can be 
interpreted as a pump-probe setup. 

(iii) Final-state electron-electron interaction in the 
continuum. After the second electron has been released 
within the short time interval T p as well, their mutual 
repulsion may redirect the electrons. 

While the dividing line between those mechanisms is 
far from sharp, the present time-dependent wave packet 
propagation allows to shed light on their relative im- 
portance since they occur on different time scales. Re- 
laxation of the ground-state correlations (i) is expected 
to occur on the time scale of the orbital period of the 
residual electron. As the remaining one-electron wave 
function will be mostly in the n = 1 and n = 2 shells, 
the timescale for this relaxation can be estimated as 
tp) ~ h/ (E 2 — Ei) w 16 as, where E n is the binding energy 
in the n-th shell of the He + ion. Therefore, ground-state 
correlations will become clearly visible only for pulses 
with durations shorter than those investigated here. The 
time scale for induced dipole polarization (ii) can be es- 
timated by the time the first electron takes to escape 



to a distance where it does not influence the remaining 
bound electron strongly. Choosing a somewhat arbitrary 
distance of 10 a. u., the time necessary for the first elec- 
tron to reach this distance after absorbing a 70 eV photon 
is about 120 as and thus of the order of the pulse lengths 
T p considered. For higher photon energies, the first elec- 
tron escapes more quickly, decreasing the importance of 
this effect. In order to verify this energy dependence, we 
have performed calculations at various photon energies 
for T p = 75 as. Fig. 4a demonstrates that for higher en- 
ergies the asymmetry of the joint angular distribution is 
indeed strongly reduced. 

Long-range Coulomb interactions in the continuum 
(iii) extend over much longer timescales which strongly 
depend on the relative emission angles and energies of the 
electrons, i.e., |ki — k 2 |. For example, for two electrons 
ejected in the same direction and with similar energies, 
the interaction will last much longer than for ejection 
in opposite directions. This can be verified by using an 
ultrashort pulse to start a two-electron wave packet in 
the continuum and observing the evolution of the joint 
angular distribution after the laser pulse is switched off 
(Fig. 4b). Directly after the pulse, the distribution of 
the electrons shows a decreased probability for ejection 
on the same side of the nucleus (primarily because of 
(ii)), but the lobes in forward and backward direction 
still mostly retain the shape expected from a dipole tran- 
sition. As time passes, continuum final-state interactions 
persist and the joint angular distribution develops a pro- 
nounced dip at equal ejection angle. The change at larger 
relative angles is almost negligible. 

One remarkable feature of the conditional angular dis- 
tribution is the persistence of the nodal plane at = 90°. 
While correlation effects strongly perturb the shape of 
the independent-particle dipolar shape, the nodal plane 
expected for the angular distribution of two electrons ab- 
sorbing one photon each is preserved. This is in contrast 
to one-photon double ionization, where only one electron 
absorbs the photon energy and electron ejection at an- 
gles normal to the polarization axis is indeed observed 
[7]. On the other hand, the conditional angular distri- 
bution in the nonsequential regime also exhibits a nodal 
plane at = 90° (cf. Fig. 3). 

Additional insights can be gained from a projection of 
the two-electron momentum onto the energy-angle plane, 

p(£i, 012, 0i=o°) = J P(E u E 2 ,n u ii 2 )dE 2 , (2) 

in coplanar geometry and for 6\ =0°. While for long 
pulses the energy of the emitted electrons is independent 
of the relative emission angle (Fig. 5c), strong energy- 
angle correlations appear for short (T p < 450 as) pulses. 
The dominant emission channel is the back-to-back emis- 
sion at equal energy sharing (E\ w 30 eV) . This corre- 
sponds precisely to the well-known Wannier ridge riding 
mode [27], previously observed in e-2e ionization pro- 
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FIG. 5: Angle-energy distribution P(Ei, #12, 0i=0°) in coplanar geometry at 70 eV photon energy for different pulse durations: 
(a) 150 as, (b) 450 as, (c) 4500 as FWHM. The side plots show the distribution integrated over, respectively, energy and angle. 



cesses [28] and also invoked in the classification of doubly- 
excited resonances [29]. Because of the large instability of 
the Wannier orbit its presence is more prevalent in break- 
up processes than in quasi-bound resonances. A second 
subdominant but equally interesting channel opens for 
short pulses at #12 = 0°, i.e., emission in the same direc- 
tion. One of the electrons is slowed down while the other 
one is accelerated. Hence, the slow electron "pushes" the 
fast electron from behind, transferring part of the energy 
absorbed from the photon field to the faster electron. 
This is the well-known post- collision interaction [30-32] 
first observed by Barker and Berry in the decay of au- 
toionizing states excited through ion impact [33]. 

In conclusion, we have shown that for attosecond 
XUV pulses the conventional scenario of "sequential" two- 
photon double ionization (TPDI) breaks down. Due to 
the small time interval between the two photoabsorption 
processes dynamical electron-electron correlations can be 
tuned by the pulse duration T p . One can view TPDI 
as an XUV-XUV pump-probe process. The angular and 
angle-energy distributions reveal the signatures of elec- 
tronic correlation induced by the Coulomb interaction 
in the intermediate bound-free complex and in the final 
state with both electrons in the continuum. For short 
pulses, two well-known scenarios, the Wannier ridge rid- 
ing mode and the post-collision interaction process, are 
simultaneously present in the two-electron emission spec- 
trum. 
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